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A simple parametric model for spherical galaxy clusters 
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ABSTRACT 

We present an analytic parametric model to describe the baryonic and dark matter distributions 
in clusters of galaxies with spherical symmetry. It is assumed that the dark matter density 
follows a Navarro, Frenk and White (NFW) profile and that the gas pressure is described by 
a generalised NFW (GNFW) profile. By further demanding hydrostatic equilibrium and that 
the local gas fraction is small throughout the cluster, one obtains unique functional forms, 
dependent on basic cluster parameters, for the radial profiles of all the properties of interest 
in the cluster. We show these profiles are consistent both with numerical simulations and 
multi-wavelength observations of clusters. We also use our model to analyse six simulated 
SZ clusters as well as A61 1 SZ data from the Arcminute Microkelvin Imager (AMI). In each 
case, we derive the radial profile of the enclosed total mass and the gas pressure and show that 
the results are in good agreement with our model prediction. 
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1 INTRODUCTION 



Analyses of observations of galaxy clusters via their X-ray emis- 
sion, gravitational lensing or Sunyaev-Zel'dovich (SZ) effect are 
often based on some parameterised cluster model for the distribu- 
tion of the cluster dark matter and the thermodynamical properties 
of its intra-cluster medium (ICM). The accuracy and robustness of 
cluster parameters derived from studies at different wavelengths de- 
pend greatly on how well the model describes the physical proper- 
ties of the cluster, and the assumptions made regarding the dynam- 
ical state of the cluster and its gas content. 

Cluster models typically assume spherical symmetry, an ideal 
gas equation-of-state, and parameterised functional forms for the 
radial distribution of two linearly-independent cluster properties, 
such as electron density and temperature (Sanderson & Ponman 
2003; Vikhlinin et al. 2005 , 2006; Laroque et al. 2006; Feroz 
et al. 2009; AMI Consortium: Zwart et al. 2011; AMI Consor- 
tium: Rodriguez-Gonzalvez et al. 2011; AMI Consortium: Hurley- 
Walker et al. 2012 and AMI Consortium: Shimwell et al. 2011 
); electron pressure and density (Nagai et al. 2007; Mroczkowski 
et al. 2009; Arnaud et al. 2010; Plagge et al. 2010 and Planck Col- 
laboration 201 Id); or electron pressure and entropy (Allison et al. 
2011; AMI Consortium: Olamaie et al. 2012). Such parameterisa- 
tions are usually supplemented by the imposition of a condition 
on the dynamical state of the cluster, most commonly hydrostatic 
equilibrium or a virial relation, together sometimes with further 
constraints, such as a constant gas fraction throughout the cluster 
and/or assorted scaling relations. All such models, with their cor- 
responding assumptions, have the potential to introduce biases in 
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the derived cluster physical parameters, depending strongly on the 
appropriateness of the assumptions made and whether or not the 
data can constrain the parameters describing the model (see AMI 
Consortium: Olamaie et al. 2012). 

A recent interesting example of a parameterised spherical 
cluster model by Mroczkowski (2011) assumes the cluster dark 
matter density to follow a parameterised Navarro, Frenk and White 
(NFW) (Navarro et al. 1997) profile and the gas pressure to be de- 
scribed by a generalised NFW (GNFW) profile (Nagai et al. 2007) 
with fixed shape parameters, both in accordance with numerical 
simulations. This model also assumes hydrostatic equilibrium and, 
crucially, a constant gas fraction (both local and enclosed) through- 
out the cluster, which is a very stringent condition. In this paper, we 
adapt this model by replacing this last condition, which is in fact in- 
consistent with the rest of the model, by the much weaker assump- 
tion that the local gas fraction throughout the cluster is small com- 
pared with unity. We show that this assumption leads to a unique 
solution for the radial dependence of all the cluster properties of 
interest (dependent on basic cluster parameters). 

Further, we analyse six simulated clusters and one real cluster 
(A611) in this frame work through their Sunyaev- Zel'dovich ef- 
fect and show that the resulting profiles agree with those predicted 
by numerical simulations and measured in multi-wavelength obser- 
vations of galaxy clusters. 



2 THE MODEL 

The first assumption in our model is a functional form for the dark 
matter density Pdm('')- Cosmological A'-body simulations suggest 
that all dark matter halos can be modelled with the spherically aver- 
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aged density profile of Navarro, Frenk and White (NFW) (Navarro 
et al. 1997) 



PdmM 



(i) 



where p s is an overall normalisation coefficient and R s is the scale 
radius where the logarithmic slope of the profile dlnp(r)/dlnr = 
-2. It is common practice also to define the halo concentration pa- 
rameter, c'200 = -if , where 7200 is the radius at which the enclosed 
mean density is 200 times the critical density at the cluster redshift. 

Our second assumption is a functional form for the gas pres- 
sure Pgas(r). Numerical simulations (Nagai et al. 2007) and X-ray 
observations of clusters of galaxies using Chandra (Vikhlinin et al. 
2005, 2006; Nagai 2006; Nagai et al. 2007 ) both show that self- 
similarity is more likely to be observed in the gas pressure profile 
than the density or temperature at large radii, i.e. up to 7-500 and 
beyond. The gas pressure is also the quantity least affected by dy- 
namical history and non-gravitational mechanisms inside the ICM. 
In particular, following Nagai et al. (2007), we assume the electron 
pressure follows the GNFW profile 
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where P ci is an overall normalisation coefficient, r p is the scale ra- 
dius. It is common to define the latter in terms of r^oo, the radius 
at which the mean enclosed density is 500 times the critical den- 
sity at the cluster redshift, and the gas concentration parameter, 
C500 = 7500 /r p . The parameters (a,b,c) describe the slopes of the 
pressure profile at r a r p , r > r p and r <K r p respectively. In the 
simplest case, we follow Arnaud et al. (2010) and fix the values of 
the gas concentration parameter and the slopes to be (C500, a, b, c) = 
(1.156, 1.0620, 5.4807, 0.3292). Arnaud et al. (2010) derived these 
values by analysing profiles of the REXCESS cluster sample ob- 
served with XMM-Newton (Bohringer et al. 2007) as well as three 
different sets of detailed numerical simulations by Borgani et al. 
(2004), Piffaretti & Valdarini (2008) and Nagai et al. (2007) that 
take into account radiative cooling, star formation, and energy feed- 
back from supernova explosions. They estimated M 50O for each 
cluster in their sample using the standard M500 - Y x scaling rela- 
tion (see Appendix B in Arnaud et al. 2010). It should be noted that 
these values are different from those used by Nagai et al. (2007), 
Mroczkowski et al. (2009), Plagge et al. (2010) and Mroczkowski 
(201 1). The Arnaud values were, however, used to analyse SZ data 
from the Planck survey data (Planck Collaboration 201 Id). 

More generally, one can allow the parameters (C500, a, b, c) to 
vary, although we will not consider this case here. Given the elec- 
tron pressure, the gas pressure is then defined by 

PUr) = -P e (r), (3) 

where /j e = 1.14m p (Mason & Myers 2000) is the mean gas mass 
per electron, p = 0.6m p is the mean mass per gas particle and m p is 
the proton mass. 

Our third assumption concerns the dynamical state of the clus- 
ter, which we take to be in hydrostatic equilibrium throughout. 
Thus, the total cluster mass internal to radius r is related to the 
gas pressure gradient at that radius by 

dfga S (r) , GM, ot (r) 

-, = -PgasM 5 • (4) 

dr r 1 
However, we note that the latest cosmological simulations of 
galaxy clusters with focus on studying the cluster outskirts (Lau 



et al. 2009; Battaglia et al. 2010; Nagai 2011; Nagai & Lau 2011; 
Parrish et al. 2012 and Battaglia et al. 2011 a,b) and observational 
studies of the clusters using the Suzaku and XMM- Newton satel- 
lites at large radii- out to the virial radius including A1795 (Bautz 
et al. 2009), PKS 0745-191 (George et al. 2009), A2204 (Reiprich 
et al. 2009), A1413 (Hoshino et al. 2010), A1689 (Kawaharada 
et al. 2010), Virgo cluster (Urban et al. 2011) and Perseus clus- 
ter (Simionescu et al. 201 1), show that the presence of random gas 
motion, gas clumping and turbulence due to the magnofhermal in- 
stability in the intracluster medium of galaxy clusters provide non- 
thermal pressure support and can introduce biases in HSE mea- 
surements of the ICM profiles and cluster mass. Hence in order to 
recover these profiles accurately we need to modify the equation of 
HSE to take into account non-thermal pressure components. How- 
ever, as the studies of this kind (to understand the physics of the 
cluster outskirts and make accurate measurements of the ICM pro- 
files in the cluster outer regions) are still ongoing, we do not study a 
modified form of HSE here. We , of course, aim to consider a more 
general form in our future analyses. 

Finally, our model is completed by assuming that the local 
gas fraction is much less than unity throughout the cluster, i.e. 

<K 1 for all r. This final assumption allows us to write 
AotW = PdmM + PgasW ~ PdmM- We emphasize that this as- 
sumption is for the gravitational part of the calculation and as we 
show in equation (6) we do not assume p gas (r) = . Thus, from (1), 
the total mass enclosed with a radius r has the analytical solution 
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Substituting this form and the expressions (2) and (3) for the gas 
pressure into the condition (4) for hydrostatic equilibrium, one may 
derive the gas density profile 
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The radial profile of the electron number density is then trivially 
obtained using n e (r) = p w (r)lp e . Assuming an ideal gas equa- 
tion of state, this in turn yields the electron temperature profile 
k B r e (r) = P c (r)/n c (r), given by 



k B r e (r) 



(4^Gp s )(^) x 
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(7) 



which is also equal to the gas temperature profile k B r gas (r). We can 
also determine the radial profile for electron entropy of the ICM. In 
the astronomy literature, for an adiabatic monatomic gas, entropy 
is defined as K s = k B T e (r)n~ 2/3 (r) which is related to the true ther- 
modynamic entropy per gas particle via 5 = f &b m (^e) + So where 
Sq is a constant (Voit 2005). 
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The only fundamental cluster property for which the radial 
profile cannot be expressed in an explicit analytical form is the gas 
mass enclosed within radius r, 



Table 1. The input parameters and ranges used in the analysis 



M gas (r) 
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For the gas density profile in (6), we have been unable to evalu- 
ate this expression analytically, and so M giis (r) must be obtained 
using numerical integration. Consequently, the enclosed gas mass 
fraction profile / eas (r) = M gils (r)/M tot (r) also cannot be writ- 
ten in closed form. It is clear, however, that the resulting f gas {r) 
will not be constant. Therefore, this contradicts the assumption of 
Mroczkowski (201 1) of / gas (r) being constant. In the next section, 
we represent the profile of / eas (r) which illustrates this point. 



3 ILLUSTRATION OF CLUSTER PROPERTIES 

In the simplest case, where a, b , c and c 500 in (2) have fixed val- 
ues, our cluster model depends only on three parameters: p s and R s 
in the NFW dark matter density profile (1) and in the pressure 
profile (2). One is, however, free to choose alternative parameters 
to define a cluster, although this choice and the priors imposed on 
the parameters can lead to very different results in the analysis of 
cluster observations (AMI Consortium: Olamaie et al. 2012). Here 
we will define clusters in terms of the parameter set M tot (r 20 o), C200, 
/gasfooo) and the redshift z, and investigate the resulting radial pro- 
files of quantities of interest in our cluster model. 

For illustration purposes, we will consider clusters at a fixed 
redshift z = 0.3. We will further assume that /gasfooo) = 0.12, 
which is reasonable since we expect the gas mass fraction to ap- 
proach the universal baryon fraction at large scales near the virial 
radius (Komatsu et al. 2011; Larson et al. 2011). We will con- 
sider a selection of 15 clusters equally spaced in the mass range 
1.0 x 10 14 M G < M tot (r 200 ) < 1.5 x 10 15 M Q . For each cluster, we 
also consider 15 values of C200 in the range 4-6, since the con- 
centration parameter shows a clear dependence on the halo mass, 
with massive halos having lower concentration parameter (Pointe- 
couteau et al. 2005; Salvador-Sole et al. 2007; Vikhlinin et al. 
2006). Throughout, we also assume a ACDM cosmology with 
n M = 0.3 , Q A = 0.7 , cr s = 0.8 , h = 0.7, w = -1 , w a = 0. 

To determine the radial profiles of quantities of interest for 
a given cluster, one must first determine the values of the model 
parameters p s , R s , r p and P ei . Since M 10[ (r 2 oo) is the total amount of 
matter internal to radius 7-200 > one may write 



An 

M tot (r 20 o) = y 4)o( 20 °Pcnt(z))- 



(9) 



Thus, for a given M tot (r 2 oo) and z, one may calculate r 20 o, and hence 
^ s = r 20o/ c 20o- The value of p s is then obtained by equating the 
input value of M tot (?- 2 oo) with the RHS of (5) evaluated at r = 7-200 > 
and is given by 



200 



Pcnt(z) 



(10) 



By equating equations (5) and (9) at r 500 , one may calculate r 500 
and hence r p = 7500 /csoo- Finally, P e j is obtained by substituting (6) 
into (8), evaluating the RHS at r = ?2oo and equating the result to 
/gastoooWtot^oo)- This yields 

P " = [jf) (G ^ )M ^ (rMo) x 



Parameter 


Value 


Aftotfooo) 


(1.0 15.0) X 10 14 M o 


C200 


(4 6) 


z 


0.3 


/gtooo) 


0.12 




Figure 1. Dark matter density profiles pdm(>')- For a given thick line, 
the thickness represents the spread in varying the halo concentration 
parameter,c 2 oo between 4 and 6. 



/"""200 






1 r' 3 dr' 







(11) 



which must be evaluated numerically. 

Tab. 1 summarises the input parameter values for our illustra- 
tive clusters. The corresponding radial profiles for various quanti- 
ties of interest are shown in Figs. 1-8. The thickness of each line 
represents the spread in halo concentration parameter correspond- 
ing to varying c 2 oo between 4 and 6-i.e. each thick line represents 
clusters with the same M lot (r 20 o) but different c 20 o- It should be 
noted that each profile is plotted out to r 200 for the corresponding 
cluster. 



4 ANALYSIS OF INTERFEROMETRIC SZ 
OBSERVATIONS 

In order to verify that our proposed model, with its corresponding 
assumptions, can describe profiles of cluster physical properties ac- 
curately, we carried out a Bayesian analysis (Feroz & Hobson 2008; 
Feroz et al. 2009a) of a set of six simulated clusters as well as A61 1 
observed through their Sunyaev-Zel'dovich (SZ) effect (Sunyaev 
& Zeldovich 1970; Birkinshaw 1999; Calrstrom, Holder & Reese 
2002) using the Arcminute Microkelvin Imager (AMI) (AMI Con- 
sortium: Zwart et al. 2008). 

The observed SZ surface brightness in the direction of electron 
reservoir may be described as 

dB v 



SU = TcMByf(v) 



8T 



(12) 



T=T CMB 



Here B v is the blackbody spectrum, T C mb = 2.73 K (Fixsen 
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Figure 8. Gas mass fraction profiles /g as (r). 

et al, 1996) is the temperature of the CMB radiation, f(v) = 
(xfr-f - 4) (1 + S(x, r e ) is the frequency dependence of thermal SZ 

signal, x = , hp is Planck's constant, v is the frequency and ks 

is Boltzmann's constant. S(x, T e ) takes into account the relativistic 
corrections in the study of the thermal SZ effect which is due to the 
presence of thermal weakly relativistic electrons in the ICM and is 
derived by solving the Kompaneets equation up to the higher orders 
(Rephaeli 1995, Itoh et al. 1998, Nozawa et al. 1998, Pointecouteau 
et al. 1998 and Challinor and Lasenby 1998). It should be noted 
that at 15 GHz (AMI observing frequency) x = 0.3 and therefore 
the relativistic correction, as shown by Rephaeli (1995), is negligi- 
ble for ^b^c < 15 keV. The dimensionless parameter y, known as 
the Comptonization parameter, is the integral of the number of col- 
lisions multiplied by the mean fractional energy change of photons 
per collision, along the line of sight 

V = —7 f n e (r)kMr)dl (13) 

where n c (r), P e (r) and T c are the electron number density, pressure 
and temperature at radius r respectively. <xx is Thomson scattering 
cross-section, m e is the electron mass, c is the speed of light and dl 
is the line element along the line of sight. It should be noted that in 
equation (13) we have used the ideal gas equation of state. 

An interferometer like AMI operating at a frequency v mea- 
sures samples from the complex visibility plane 7 v (u). These are 
given by a weighted Fourier transform of the surface brightness 
/„(x), namely 

7 F (u) = J A y (x)/,,(x) exp(27riu • x)dx, (15) 

where x is the position relative to the phase centre, A v (x) is the 
(power) primary beam of the antennas at observing frequency v 
(normalised to unity at its peak) and u is the baseline vector in 
units of wavelength. 

Further details of our Bayesian methodology, modelling in- 
terferometric SZ data, primordial CMB anisotropies, and resolved 
and unresolved radio point-source models are given in Hobson & 
Maisinger (2002), Feroz & Hobson (2008) and (2009 a,b), AMI 
Consortium: Davies et al. (2011) and AMI Consortium: Olamaie 
et al. (2012). 



A parametric model for galaxy clusters 5 

Table 2. The input parameters used to generate simulated clusters . 



Cluster 


M lot (r 20 o)10 14 M 


C200 




/gfooo) 


clsiml 


5.0 


4.0 


0.3 


0.12 


clsim2 


5.0 


4.3 


0.3 


0.12 


clsim3 


5.0 


4.5 


0.3 


0.12 


clsim4 


5.0 


5.0 


0.3 


0.12 


clsim5 


5.0 


5.5 


0.3 


0.12 


clsim6 


5.0 


6.0 


0.3 


0.12 



In generating simulated SZ skies and observing them with a 
model AMI SA, we have used the methods outlined in Hobson & 
Maisinger (2002), Grainge et al. (2002), Feroz et al. (2009b) and 
AMI Consortium: Olamaie et al. (2012). 

Generating a simulated cluster SZ signal using the model de- 
scribed in Sections 2 and 3 requires the input parameters M lot (r2oo), 
C200* z and / s , as (''20o) listed in Tab. 1; this set of parameters fully 
describes the Comptonization y parameter. Tab. 2 summarises the 
input parameters used to generate six simulated SZ clusters. The 
simulated clusters all have the same M to ,(f2oo), z and / gas (r2oo) and 
the only parameter that varies from cluster to cluster is the halo 
concentration parameter, C2oo- 

A61 1 is a rich cluster at redshift z = 0.288 and has been stud- 
ied through its X-ray emission, strong lensing, weak lensing and SZ 
effect (Schmidt & Allen 2007; Romano et al. 2010; Donnarumma 
et al. 201 1 and AMI Consortium: Shimwell et al. 201 l).These stud- 
ies suggest that there is no significant contamination from radio 
sources and there is no evidence for a radio halo associated with 
A611 (Venturi et al. 2008). The SZ signal (decrement) on the AMI 
map appears circular, (fig. 2 in AMI Consortium: Shimwell et al. 
2011) in agreement with the X-ray surface brightness from the 
Chandra archive data (fig. 2 in AMI Consortium: Shimwell et al. 
2011 and fig. 1 in Donnarumma et al. 2011), which also appears 
to be smooth and whose peak coincides with the position of the 
brightest cluster galaxy and the SZ peak. These results are a strong 
indication that the cluster is relaxed. Moreover, the absence of radio 
halos in the cluster which are major sources of the presence of non- 
thermal mechanisms in the galaxy clusters ( Brunetti et al. 2009) 
makes A611 an ideal cluster candidate for our analysis as it satis- 
fies both assumptions of spherical symmetry and thermal pressure 
support in equation of the HSE. 

Details of AMI pointed observation towards the cluster, data 
reduction pipeline and mapping are described in AMI Consortium: 
Shimwell et al. (2011) and in here we focus on the Bayesian anal- 
ysis of the clusters using our model. 

The sampling parameters in our Bayesian analysis are C = 
(-tc,yc,c 2 oo,Mot(''20o),/g(''20o),z), where x c and y c are cluster pro- 
jected position on the sky. We further assume that the priors on 
sampling parameters are separable (Feroz et al. 2009b) such that 

7r(0 t .) = 7r(x c )^(y c )^(c2oo)?r(WT('-20o))?r(/g('-20o))7r(z)- (16) 

We use Gaussian priors on cluster position parameters, centred on 
the pointing centre and with standard-deviation of 1 arcmin. We 
adopt uniform priors on C200 and a S function prior on redshift z. 
The prior on M lot (r 2 oo) is taken to be uniform in logM in the range 
M min = 10 14 M Q to M max = 6 x 10 15 M and the prior of / gas (r 2 oo) is 
set to be a Gaussian centred at the / gas = 0. 12 with a width of 0.02. 
A summary of the priors and their ranges are presented in Tab. 3. 

In order to understand the underlying biases and constraints 
imposed by the priors and the model assumptions, we first study 
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Table 3. Summary of the priors on the sampling parameters. Note that 
N(ji, <t) represents a Gaussian probability distribution with mean /i and stan- 
dard deviation of <x and U(a, b) represents a uniform distribution between a 
and b. 



Parameter 


Prior 


x c > yc 


N(0 , 60)" 


t'200 


U(l , 10) 


log M tot (r 2 oo) 


C/(14 , 15.8) M G 


/gasteoo) 


W(0.12 , 0.02) 




-200 200 -200 200 5 10 2 4 6 0.1 0.2 



x /arcsec y /arcsec c 200 i g Jr 2m ) 

Figure 9. 1 -D and 2-D marginalised posterior distributions of sampling pa- 
rameters with no data . 

our methodology in the absence of data. This can be carried out by 
setting the likelihood to a constant value and hence the algorithm 
just explores the prior space. Along with the analysis done using the 
simulated AMI data, this approach reveals the constraints that mea- 
surements of the SZ signal place on the cluster physical parameters 
and the robustness of the assumptions made. Fig. 9 represents 1-D 
and 2-D marginalised posterior distributions of a prior-only anal- 
ysis for each of the sampling parameters in our model. The plots 
show that we correctly recover the assumed prior probability dis- 
tributions of the sampling parameters in the absence of SZ data. 

Fig. 10 shows 1-D and 2-D marginalised posterior distribu- 
tions of sampling parameters for the first simulated SZ cluster 
data, with vertical lines representing the true parameter values and 
Fig. 1 1 shows the results of the analysis for A61 1. From the plots 
we notice that the model, along with its corresponding assumptions, 
can constrain cluster position and M tol (r 2 oo). but C200 remains rela- 
tively unconstrained. We also notice the weak negative degeneracy 
between M lot (r 20 o) an d c 2 oo as we expect in high mass halos, be- 
tween l.Ox 10 14 M Q and 15.0 x 10 14 M Q (Pointecouteau et al. 2005; 
Salvador-Sole rt al. 2007; Rudd, Zentner & Kravtsov 2008; Bhat- 
tacharya, Habib, & Heitmann 2011). From our analysis we find 
M tol (r 20 o) = (5.3 ± 2.6) x 10 14 M o and r 20 o = (1.5 ± 0.2) Mpc for 
the simulated cluster and M lot (r 20 o) = (8.6 ± 1.4) x 10 14 M o and 
r 20Q = (1.7±0.1)Mpcfor A611. 

Figs. 12 and 13 present 1-D marginalised posterior distribu- 




Figure 10. 1-D and 2-D marginalised posterior distributions of sampling 
parameters for the first simulated cluster . 




Figure 11. 1-D and 2-D marginalised posterior distributions of sampling 
parameters for A61 1 . 



tions of the model parameters (i.e. p s , R s , r p and P e i) f° r the first 
simulated cluster and A61 1 respectively. We note that although our 
data can constrain P ei and r p , p s and R s are not well constrained as 
they depend strongly on the relatively unconstrained cluster con- 
centration parameter c 2 oo- We use the best-fit values of these four 
parameters given in Tabs. 4 and 5 to determine the radial profiles 
of the clusters physical parameters according to our model. 

As the SZ surface brightness is proportional to the line-of- 
sight integral of the electron pressure, (equations 12 and 13) SZ 
analysis of galaxy clusters provides a direct measurement of the 
pressure distribution of the ICM. Moreover, the integral of the 
Comptonization y parameter over the solid angle CI subtended by 




.2 



Figure 12. 1-D marginalised posterior distributions of the model parameters 
for the first simulated cluster. Green vertical lines are the derived values of 
the model parameters for a cluster defined by the input parameters given in 
the first row of Tab. 2 



Table 4. Simulated cluster model parameters estimated (mean and standard 
deviation) 



Parameter 


M 


(T 


p s (M Mpc- 3 ) 


2.52 X 10 15 


2.31 x 10 15 


R s (kpc) 


402.79 


306.83 


fp(kpc) 


846.93 


95.88 


PeiCkeVnr 3 ) 


1.54 x 10 4 


3.08 x 10 3 



the cluster (Yg Z ) is proportional to the volume integral of the gas 
pressure. It is thus a good estimate for the total thermal energy con- 
tent of the cluster and its mass (see e.g. Bartlett & Silk 1994). Y$z 
parameter in both cylindrical and spherical geometries may be de- 
scribed as 

Y<~l(R) = f dl f P c (r)27T5ds (17) 

Y spb (r) = """"T f ^('"')4^' 2 dr' (18) 
m e c i Jo 

where R is the projected radius of the cluster on the sky. In this con- 
Table S. Best-fit values of model parameters estimated (mean and standard 
deviation) for A61 1. 



Parameter 


/' 


IT 


p s (M Mpc- 3 ) 


2.87 x 10 15 


2.07 x 10 15 


R s (kpc) 


361.80 


209.31 


r p (kpc) 


1021.27 


55.57 


P ei (keVnr 3 ) 


2.41 x 10 4 


2.59 x 10 3 
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Figure 13. 1-D marginalised posterior distributions of the model parameters 
for A611. 




0.5 1 1.5 2 2.5 



r (Mpc) 

Figure 14. Profile of M tot versus r for simulated clusters. 

text we determined the radial profiles of M tol and P c as a function 
of r for six simulated clusters and A61 1 (figs. 14-17). In all figures, 
the background thick line shows the model prediction of the pro- 
files of the clusters with the same mass as the clusters been analysed 
but vary in c 2 oo as was illustrated in figs. 1-8. We have plotted the 
radial profiles of the corresponding cluster properties with coloured 
* and « in case of simulated clusters and black * for A611. From 
the plots it is obvious that the radial trend of the clusters profiles 
are all consistent with our model prediction. 

5 DISCUSSION AND CONCLUSIONS 

We have studied the combination of NFW dark matter and GNFW 
gas pressure profiles within the hierarchical structure formation 



scenario (Kaiser 1986) to derive the radial distribution of the clus- 
ter properties, assuming spherical symmetry, hydrostatic equilib- 
rium, and that the local gas fraction throughout is small compared 
to unity. 

Figs. 1 and 2 represent the dark matter density and enclosed 
mass profiles (the latter approximating the total enclosed mass pro- 
file under our assumption that <sc 1 throughout). These results 
are based on the assumption of an NFW dark matter density profile, 
arising from the notion that the largest virialised structures form via 
gravitational collapse and merging. The profile has proved to be a 
good fit to the relaxed dark matter halos in high resolution N-body 
simulations down to 1 % of the virial radius, and optical and X-ray 
observations of galaxy clusters both indicate that the profile is a 
good representation of the underlying cluster mass profile outside 
the core (Carlberg et al. 1997; Pratt & Arnaud 2002). 

The electron pressure profiles are shown in Fig. 3, which are 
assumed to follow a GNFW profile. They exhibit self-similarity at 
the larger radii as they approach r 2 oo and show dependency on the 
cluster mass. These behaviours are expected and resemble the pro- 
files observed in a wide range of clusters (Holder et al, 2007; Nagai 
et al. 2007; Arnaud et al. 2010; Plagge et al.2010; Mroczkowski 
et al. 2009) indicating that the pressure is least affected by non- 
gravitational phenomena in the ICM. This is important, in partic- 
ular, in the analysis of Sunyaev-Zel'dovich observations of clus- 
ters, which essentially measure the line-of-sight integral of the ICM 
pressure through the cluster. 

The derived gas density profiles p gas (r) = jj. e ii e (r) are shown in 
Fig. 4, and reproduce all the main features observed both in numer- 
ical simulations and in real clusters (Sanderson et al. 2003; Bor- 
gani et al. 2004; Vikhlinin et al. 2005, 2006; Nagai et al. 2007; 
McCarthy et al. 2008). In particular, the profiles exhibit steepening 
at large radii, a power-law cusp at small radii (resulting from the 
fact that gas cooling and star formation processes have been taken 
into account in deriving the GNFW pressure profile) and a change 
of slope at intermediate radii. We also note that lower temperature 
clusters have lower gas density. The derived analytical expression 
for the gas density can thus model both the inner and outer regions 
of the clusters. 

The derived electron/gas temperature profiles of the clusters 
are shown in Fig. 5. All of them have similar positive slopes up to 
r ~ 0. 1 Mpc for the most massive clusters and have a broad peak 
around this region. The significant drop in the temperature in the 
innermost region is again because of taking into account the pres- 
ence of radiative cooling mechanisms in deriving the GNFW pres- 
sure profile (Borgani et al. 2004; Vikhlinin et al. 2005, 2006; Pratt 
et al. 2007; Leccardi & Molendi 2008). In particular, we note that 
the clusters do not have isothermal cores. It should be pointed out, 
however, that real cluster data and current high-resolution simula- 
tions display complex temperature structures, which are the result 
of merging subgroups or supersonic accretion which heats the gas 
across the shock front where the assumption of hydrostatic equilib- 
rium clearly breaks. Nonetheless, our derived temperature profile 
describes the general features of the ICM well, within our assump- 
tions. 

Fig. 6 shows the entropy profiles, which clearly show that the 
entropy depends on temperature and therefore the cluster mass. 
Moreover, the entropy profiles approach self-similarity as the ra- 
dius approaches r 2 oo> showing a scaling power-law distribution 
which is predicted in the models based on spherical gas accre- 
tion within a NFW dark matter halo (Tozzi & Norman 2001). This 
demonstrates that gravity dominates the ICM thermodynamics in 
the outer regions of clusters. These behaviours have already been 
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noted in the cluster numerical simulations (Kay et al. 2004; Bor- 
gani et al. 2004; Voit, Kay & Bryan 2005) and have also been ob- 
served in a large sample of galaxy clusters (Ponman, Sanderson and 
Finoguenov 2003; Piffaretti et al. 2005; McCarthy et al. 2008; Pratt 
et al. 2010). In the inner regions, on the other hand, the entropy 
profiles are clearly affected by the non-gravitational processes that 
have been considered in deriving the GNFW pressure profile. 

Fig. 7 presents the enclosed gas mass profiles which increas- 
ing with radius but with different slopes and fig. (8) shows the de- 
rived gas mass fraction profiles, which also exhibit a significant 
increase with radius, hence implying that / gas can not be constant 
throughout the cluster as assumed by Mroczkowski (201 1). Indeed, 
such an assumption is inconsistent with our other model assump- 
tions as they lead to / gas being a function of r. The profiles also 
show a pronounced dependency on the cluster mass, reflecting the 
dependency on the temperature as expected both from numerical 
simulations and X-ray observations of galaxy clusters using XMM 
and Chandra satellites (Ettori et al. 2004; Allen et al. 2004; Sa- 
dat et al. 2005; Vikhlinin et al. 2005, 2006; LaRoque et al. 2006; 
McCarthy et al. 2007; Afshordi et al.2007). 

Moreover, by numerically exploring the probability distribu- 
tions of the cluster parameters given simulated interferometric SZ 
data in the context of Bayesian methods, and assuming our model 
with its corresponding assumptions, we investigate the capability of 
this model and analysis to return the simulated cluster input quan- 
tities. We find that simulated cluster physical parameters are well- 
constrained except C200 which is relatively unconstrained. We can 
also recover the true values of the simulated clusters. In particu- 
lar, the mean cluster total mass estimate M tot (r 20 o) and r 2 m for the 
first simulated cluster are: M tol (r 2 „o) = (5.1 ± 1.7) x 10 14 M o and 
'mo = (1-5 ± 0.2) Mpc and the corresponding true values of the 
simulated cluster are: M tol (r 2 oo) = 5 x 10 14 M Q and r 2 oo =1-5 Mpc. 
We determine the best-fit values of the parameters describing our 
model, i.e. p s , R s , r p and P ci , and hence calculate profiles of cluster 
total mass and gas pressure as determined using SZ data. We show 
that these profiles are consistent with our model. 

We then repeat the analysis for a real cluster (A611) ob- 
served through its SZ effect with AMI. For this cluster, We find 
MUriao) = (8.6 ± 1.4) x 10 14 M o and r 200 = (1.7 ± 0.1) Mpc. A611 
has previously been studied in different wave-bands. For example, 
Schmidt & Allen (2007) analysed Chandra data of A61 1 and found 
Moi(''20o) ~ 8 x 10 14 M o and r 2 m = 1.7 Mpc. Donnarumma et al. 
(201 1) also studied Chandra X-ray data of A61 1 with different as- 
sumptions on background and metallicity. Their estimates of the 
cluster total mass vary from M tot (r 20 o) = (9.32 ± 1.39) x 10 14 M o 
for r 20 o ~ 1.8Mpc to M tol (r 200 ) = (11.1 1 ± 2.06) x 10 14 M o for 
r 2 oo ~ 1.96 Mpc. They also carried out a strong lensing analysis 
of the cluster and found the mass estimates vary from M tot (r 200 ) = 
(4.68±0.31)xl0 14 M G forr 200 « 1.4 Mpc to M tot (r 200 ) = 6.32+°;^ x 
10 14 M for r 2 oo ~ 1.5 Mpc when using different techniques. From 
weak lensing study of the cluster, Romano et al. (2010) find that the 
cluster total mass within radius of 1.5Mpc is (8 + 3) x 10 14 M Q from 
the aperture mass technique and (5 ± 1) x 10 14 M assuming para- 
metric models. Our previous SZ analysis of A611 using isother- 
mal p model (AMI Consortium: Shimwell et al. 2011) resulted in 
M tot (r 20 o) = (5.7 ± 1.1) x 10 14 M Q and r 200 = (1.6 ±0.1) Mpc. Com- 
paring the results of these studies with our analysis reveals that our 
results are in good agreement with the results of X-ray and weak 
lensing analyses of A611 but strong lensing and our previous SZ 
analyses of the cluster find a lower cluster mass which might be 
due to the extrapolating the strong lensing results in the outer spa- 
tial range as has been pointed out by Donnarumma et al. (2011) 



and also the assumption of isothermality in our previous SZ study 
of the cluster. 

We conclude that our proposed simple model for spherical 
galaxy clusters leads to realistic radial profiles for all the prop- 
erties of interest, and hence may prove useful in the analysis of 
multi-wavelength cluster observations. An obvious future avenue 
for research, which we will explore in a follow-up paper to this 
letter, is to iterate the solution we have obtained by inserting the 
derived p gas (r) in (6) back into the expression for the total density 
PtoiM = PdmW + PgasM* recalculating the form of the other vari- 
ables and repeating this process until convergence is established. In 
so doing, one might hope to obtain an even more realistic cluster 
model, but at the cost of losing a simple analytical formulation. 
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